Whole-Genome Sequencing Analysis of Human Metabolome in Multi-Ethnic Populations.

Circulating metabolite levels may reflect the state of the human organism in health and disease, however, the genetic architecture of metabolites is not fully understood. We have performed a whole-genome sequencing association analysis of both common and rare variants in up to 11,840 multi-ethnic participants from five studies with up to 1666 circulating metabolites. We have discovered 1985 novel variant-metabolite associations, and validated 761 locus-metabolite associations reported previously. Seventy-nine novel variant-metabolite associations have been replicated, including three genetic loci located on the X chromosome that have demonstrated its involvement in metabolic regulation. Gene-based analysis have provided further support for seven metabolite-replicated loci pairs and their biologically plausible genes. Among those novel replicated variant-metabolite pairs, follow-up analyses have revealed that 26 metabolites have colocalized with 21 tissues, seven metabolite-disease outcome associations have been putatively causal, and 7 metabolites might be regulated by plasma protein levels. Our results have depicted the genetic contribution to circulating metabolite levels, providing additional insights into understanding human disease.

A multi-ethnic polygenic risk score is associated with hypertension prevalence and progression throughout adulthood.

In a multi-stage analysis of 52,436 individuals aged 17-90 across diverse cohorts and biobanks, we train, test, and evaluate a polygenic risk score (PRS) for hypertension risk and progression. The PRS is trained using genome-wide association studies (GWAS) for systolic, diastolic blood pressure, and hypertension, respectively. For each trait, PRS is selected by optimizing the coefficient of variation (CV) across estimated effect sizes from multiple potential PRS using the same GWAS, after which the 3 trait-specific PRSs are combined via an unweighted sum called "PRSsum", forming the HTN-PRS. The HTN-PRS is associated with both prevalent and incident hypertension at 4-6 years of follow up. This association is further confirmed in age-stratified analysis. In an independent biobank of 40,201 individuals, the HTN-PRS is confirmed to be predictive of increased risk for coronary artery disease, ischemic stroke, type 2 diabetes, and chronic kidney disease.

Host and gut microbial tryptophan metabolism and type 2 diabetes: an integrative analysis of host genetics, diet, gut microbiome and circulating metabolites in cohort studies.

OBJECTIVE: Tryptophan can be catabolised to various metabolites through host kynurenine and microbial indole pathways. We aimed to examine relationships of host and microbial tryptophan metabolites with incident type 2 diabetes (T2D), host genetics, diet and gut microbiota. METHOD: We analysed associations between circulating levels of 11 tryptophan metabolites and incident T2D in 9180 participants of diverse racial/ethnic backgrounds from five cohorts. We examined host genome-wide variants, dietary intake and gut microbiome associated with these metabolites. RESULTS: Tryptophan, four kynurenine-pathway metabolites (kynurenine, kynurenate, xanthurenate and quinolinate) and indolelactate were positively associated with T2D risk, while indolepropionate was inversely associated with T2D risk. We identified multiple host genetic variants, dietary factors, gut bacteria and their potential interplay associated with these T2D-relaetd metabolites. Intakes of fibre-rich foods, but not protein/tryptophan-rich foods, were the dietary factors most strongly associated with tryptophan metabolites. The fibre-indolepropionate association was partially explained by indolepropionate-associated gut bacteria, mostly fibre-using Firmicutes. We identified a novel association between a host functional LCT variant (determining lactase persistence) and serum indolepropionate, which might be related to a host gene-diet interaction on gut Bifidobacterium, a probiotic bacterium significantly associated with indolepropionate independent of other fibre-related bacteria. Higher milk intake was associated with higher levels of gut Bifidobacterium and serum indolepropionate only among genetically lactase non-persistent individuals. CONCLUSION: Higher milk intake among lactase non-persistent individuals, and higher fibre intake were associated with a favourable profile of circulating tryptophan metabolites for T2D, potentially through the host-microbial cross-talk shifting tryptophan metabolism toward gut microbial indolepropionate production.

Exome sequence association study of levels and longitudinal change of cardiovascular risk factor phenotypes in European Americans and African Americans from the Atherosclerosis Risk in Communities Study.

Cardiovascular disease (CVD) is responsible for 31% of all deaths worldwide. Among CVD risk factors are age, race, increased systolic blood pressure (BP), and dyslipidemia. Both BP and blood lipids levels change with age, with a dose-dependent relationship between the cumulative exposure to hyperlipidemia and the risk of CVD. We performed an exome sequence association study using longitudinal data with up to 7805 European Americans (EAs) and 3171 African Americans (AAs) from the Atherosclerosis Risk in Communities (ARIC) study. We assessed associations of common (minor allele frequency > 5%) nonsynonymous and splice-site variants and gene-based sets of rare variants with levels and with longitudinal change of seven CVD risk factor phenotypes (BP traits: systolic BP, diastolic BP, pulse pressure; lipids traits: triglycerides, total cholesterol, high-density lipoprotein cholesterol [HDL-C], low-density lipoprotein cholesterol [LDL-C]). Furthermore, we investigated the relationship of the identified variants and genes with select CVD endpoints. We identified two novel genes: DCLK3 associated with the change of HDL-C levels in AAs and RAB7L1 associated with the change of LDL-C levels in EAs. RAB7L1 is further associated with an increased risk of heart failure in ARIC EAs. Investigation of the contribution of genetic factors to the longitudinal change of CVD risk factor phenotypes promotes our understanding of the etiology of CVD outcomes, stressing the importance of incorporating the longitudinal structure of the cohort data in future analyses.

Genome-wide association study of serum metabolites in the African American Study of Kidney Disease and Hypertension.

The genome-wide association study (GWAS) is a powerful means to study genetic determinants of disease traits and generate insights into disease pathophysiology. To date, few GWAS of circulating metabolite levels have been performed in African Americans with chronic kidney disease. Hypothesizing that novel genetic-metabolite associations may be identified in a unique population of African Americans with a lower glomerular filtration rate (GFR), we conducted a GWAS of 652 serum metabolites in 619 participants (mean measured glomerular filtration rate 45 mL/min/1.73m2) in the African American Study of Kidney Disease and Hypertension, a clinical trial of blood pressure lowering and antihypertensive medication in African Americans with chronic kidney disease. We identified 42 significant variant metabolite associations. Twenty associations had been previously identified in published GWAS, and eleven novel associations were replicated in a separate cohort of 818 African Americans with genetic and metabolomic data from the Atherosclerosis Risk in Communities Study. The replicated novel variant-metabolite associations comprised eight metabolites and eleven distinct genomic loci. Nine of the replicated associations represented clear enzyme-metabolite interactions, with high expression in the kidneys as well as the liver. Three loci (ACY1, ACY3, and NAT8) were associated with a common pool of metabolites, acetylated amino acids, but with different individual affinities. Thus, extensive metabolite profiling in an African American population with chronic kidney disease aided identification of novel genome-wide metabolite associations, providing clues about substrate specificity and the key roles of enzymes in modulating systemic levels of metabolites.

Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals.

Genetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to ~1.3 million participants, we discovered 106 new BP-associated genomic regions and 87 rare (minor allele frequency ≤ 0.01) variant BP associations (P < 5 × 10-8), of which 32 were in new BP-associated loci and 55 were independent BP-associated single-nucleotide variants within known BP-associated regions. Average effects of rare variants (44% coding) were ~8 times larger than common variant effects and indicate potential candidate causal genes at new and known loci (for example, GATA5 and PLCB3). BP-associated variants (including rare and common) were enriched in regions of active chromatin in fetal tissues, potentially linking fetal development with BP regulation in later life. Multivariable Mendelian randomization suggested possible inverse effects of elevated systolic and diastolic BP on large artery stroke. Our study demonstrates the utility of rare-variant analyses for identifying candidate genes and the results highlight potential therapeutic targets.

A Genome-wide Association Study Discovers 46 Loci of the Human Metabolome in the Hispanic Community Health Study/Study of Latinos.

Variation in levels of the human metabolome reflect changes in homeostasis, providing a window into health and disease. The genetic impact on circulating metabolites in Hispanics, a population with high cardiometabolic disease burden, is largely unknown. We conducted genome-wide association analyses on 640 circulating metabolites in 3,926 Hispanic Community Health Study/Study of Latinos participants. The estimated heritability for 640 metabolites ranged between 0%-54% with a median at 2.5%. We discovered 46 variant-metabolite pairs (p value < 1.2 × 10-10, minor allele frequency ≥ 1%, proportion of variance explained [PEV] mean = 3.4%, PEVrange = 1%-22%) with generalized effects in two population-based studies and confirmed 301 known locus-metabolite associations. Half of the identified variants with generalized effect were located in genes, including five nonsynonymous variants. We identified co-localization with the expression quantitative trait loci at 105 discovered and 151 known loci-metabolites sets. rs5855544, upstream of SLC51A, was associated with higher levels of three steroid sulfates and co-localized with expression levels of SLC51A in several tissues. Mendelian randomization (MR) analysis identified several metabolites associated with coronary heart disease (CHD) and type 2 diabetes. For example, two variants located in or near CYP4F2 (rs2108622 and rs79400241, respectively), involved in vitamin E metabolism, were associated with the levels of octadecanedioate and vitamin E metabolites (gamma-CEHC and gamma-CEHC glucuronide); MR analysis showed that genetically high levels of these metabolites were associated with lower odds of CHD. Our findings document the genetic architecture of circulating metabolites in an underrepresented Hispanic/Latino community, shedding light on disease etiology.

Unraveling the functional role of the orphan solute carrier, SLC22A24 in the transport of steroid conjugates through metabolomic and genome-wide association studies.

Variation in steroid hormone levels has wide implications for health and disease. The genes encoding the proteins involved in steroid disposition represent key determinants of interindividual variation in steroid levels and ultimately, their effects. Beginning with metabolomic data from genome-wide association studies (GWAS), we observed that genetic variants in the orphan transporter, SLC22A24 were significantly associated with levels of androsterone glucuronide and etiocholanolone glucuronide (sentinel SNPs p-value <1x10-30). In cells over-expressing human or various mammalian orthologs of SLC22A24, we showed that steroid conjugates and bile acids were substrates of the transporter. Phylogenetic, genomic, and transcriptomic analyses suggested that SLC22A24 has a specialized role in the kidney and appears to function in the reabsorption of organic anions, and in particular, anionic steroids. Phenome-wide analysis showed that functional variants of SLC22A24 are associated with human disease such as cardiovascular diseases and acne, which have been linked to dysregulated steroid metabolism. Collectively, these functional genomic studies reveal a previously uncharacterized protein involved in steroid homeostasis, opening up new possibilities for SLC22A24 as a pharmacological target for regulating steroid levels.

Remnant-Like Particle Cholesterol, Low-Density Lipoprotein Triglycerides, and Incident Cardiovascular Disease.

BACKGROUND: Hypertriglyceridemia is associated with increased remnant-like particle cholesterol (RLP-C) and triglycerides in low-density lipoprotein (LDL-TG). Recent studies have focused on atherogenicity of RLP-C, with few data on LDL-TG. OBJECTIVES: The aim of this study was to examine associations of RLP-C and LDL-TG with incident cardiovascular disease (CVD) events and genetic variants in the ARIC (Atherosclerosis Risk In Communities) study. METHODS: Fasting plasma RLP-C and LDL-TG levels were measured in 9,334 men and women without prevalent CVD. Participants were followed for incident CVD events (coronary heart disease and ischemic stroke) for up to 16 years. Associations between LDL-TG and RLP-C levels and genetic variants were assessed by whole-exome sequencing using single-variant analysis for common variants and gene-based burden tests for rare variants; both an unbiased and a candidate gene approach were explored. RESULTS: RLP-C and LDL-TG levels were correlated with triglyceride levels (r = 0.85 and r = 0.64, p < 0.0001). In minimally adjusted analyses, RLP-C and LDL-TG were associated with CVD risk, but in models adjusted for traditional risk factors including lipids, only LDL-TG was associated with incident CHD (hazard ratio: 1.28; 95% confidence interval: 1.10 to 1.50) and stroke (hazard ratio: 1.47; 95% confidence interval: 1.13 to 1.92). A common APOE variant, rs7412, had the strongest association with LDL-TG and RLP-C (p < 5 × 10-8). CONCLUSIONS: RLP-C and LDL-TG levels were predictive of CVD and associated with APOE variants. LDL-TG may represent a marker of dysfunctional remnant lipoprotein metabolism associated with increased CVD risk. Further research is needed to determine whether LDL-TG plays a causal role in CVD and may be a target for therapy.

Sequence-Based Analysis of Lipid-Related Metabolites in a Multiethnic Study.

Small molecule lipid-related metabolites are important components of fatty acid and steroid metabolism-two important contributors to human health. This study investigated the extent to which rare and common genetic variants spanning the human genome influence the lipid-related metabolome. Sequence data from 1552 European-Americans (EA) and 1872 African-Americans (AA) were analyzed to examine the impact of common and rare variants on the levels of 102 circulating lipid-related metabolites measured by a combination of chromatography and mass spectroscopy. We conducted single variant tests [minor allele frequency (MAF) > 5%, statistical significance P-value ≤ 2.45 × 10-10] and tests aggregating rare variants (MAF ≤ 5%) across multiple genomic motifs, such as coding regions and regulatory domains, and sliding windows. Multiethnic meta-analyses detected 53 lipid-related metabolites-locus pairs, which were inspected for evidence of consistent signal between the two ethnic groups. Thirty-eight lipid-related metabolite-genomic region associations were consistent across ethnicities, among which seven were novel. The regions contain genes that are related to metabolite transport (SLC10A1) and metabolism (SCD, FDX1, UGT2B15, and FADS2). Six of the seven novel findings lie in expression quantitative trait loci affecting the expression levels of 14 surrounding genes in multiple tissues. Imputed expression levels of 10 of the affected genes were associated with four corresponding lipid-related traits in at least one tissue. Our findings offer valuable insight into circulating lipid-related metabolite regulation in a multiethnic population.

Correction: DeltaNp63alpha-Mediated Induction of Epidermal Growth Factor Receptor Promotes Pancreatic Cancer Cell Growth and Chemoresistance.

[This corrects the article DOI: 10.1371/journal.pone.0026815.].

Whole-genome sequencing study of serum peptide levels: the Atherosclerosis Risk in Communities study.

Oligopeptides are important markers of protein metabolism, as they are cleaved from larger polypeptides and proteins. Genetic association studies may help elucidate their origin and function. In 1,552 European Americans and 1,872 African Americans of the Atherosclerosis Risk in Communities study, we performed whole-genome and whole-exome sequencing and measured serum levels of 25 peptides. Common variants (minor allele frequency > 5%) were analysed individually. We grouped low-frequency variants (minor allele frequency ≤ 5%) by a genome-wide sliding window using region-based aggregate tests. Furthermore, low-frequency regulatory variants were grouped by gene, as were functional coding variants. All analyses were performed separately in each ancestry group and then meta-analysed. We identified 22 common variant associations with peptide levels (P-value < 4.2 × 10-10), including 16 novel gene-peptide pairs. Notably, variants in kinin-kallikrein genes KNG1, F12, KLKB1, and ACE were associated with several different peptides. Variants in KLKB1 and ACE were associated with a fragment of complement component 3f. Both common variants and low-frequency coding variants in CPN1 were associated with a fibrinogen cleavage peptide. Four sliding windows were significantly associated with peptide levels (P-value < 4.2 × 10-10). Our results highlight the importance of the kinin-kallikrein system in the regulation of serum peptide levels, strengthen the evidence for a broad link between the kinin-kallikrein and complement systems, and suggest a role of CPN1 in the conversion of fibrinogen to fibrin.

Whole genome sequence analysis of serum amino acid levels.

BACKGROUND: Blood levels of amino acids are important biomarkers of disease and are influenced by synthesis, protein degradation, and gene-environment interactions. Whole genome sequence analysis of amino acid levels may establish a paradigm for analyzing quantitative risk factors. RESULTS: In a discovery cohort of 1872 African Americans and a replication cohort of 1552 European Americans we sequenced exons and whole genomes and measured serum levels of 70 amino acids. Rare and low-frequency variants (minor allele frequency ≤5%) were analyzed by three types of aggregating motifs defined by gene exons, regulatory regions, or genome-wide sliding windows. Common variants (minor allele frequency >5%) were analyzed individually. Over all four analysis strategies, 14 gene-amino acid associations were identified and replicated. The 14 loci accounted for an average of 1.8% of the variance in amino acid levels, which ranged from 0.4 to 9.7%. Among the identified locus-amino acid pairs, four are novel and six have been reported to underlie known Mendelian conditions. These results suggest that there may be substantial genetic effects on amino acid levels in the general population that may underlie inborn errors of metabolism. We also identify a predicted promoter variant in AGA (the gene that encodes aspartylglucosaminidase) that is significantly associated with asparagine levels, with an effect that is independent of any observed coding variants. CONCLUSIONS: These data provide insights into genetic influences on circulating amino acid levels by integrating -omic technologies in a multi-ethnic population. The results also help establish a paradigm for whole genome sequence analysis of quantitative traits.

DeltaNp63alpha-mediated induction of epidermal growth factor receptor promotes pancreatic cancer cell growth and chemoresistance.

Pancreatic ductal adenocarcinoma (PDAC) is highly resistant to current chemotherapy regimens, in part due to alterations in the p53 tumor suppressor pathway. p53 homolog p63 is a transcription factor essential for the development and differentiation of epithelial surfaces. However its function in cancer is controversial and its role in PDAC is not known. We discovered that ΔNp63α was the predominantly expressed p63 variant in pancreatic cancer cell lines. ΔNp63α protein and mRNA levels were high in T3M4, BxPC3 and COLO-357 pancreatic cancer cells and low in ASPC-1 and PANC-1 cells. Overexpression of ΔNp63α in PANC-1 cells and shRNA-mediated knockdown in T3M4 cells indicated that ΔNp63α promoted anchorage-dependent and -independent growth, motility and invasion, and enhanced resistance to cisplatin-induced apoptosis. Epidermal growth factor receptor (EGFR) signaling pathways contribute to the biological aggressiveness of PDAC, and we found that the motogenic effects of ΔNp63α were augmented in presence of EGF. Ectopic expression of ΔNp63α resulted in upregulation of EGFR and ß1-integrin in PANC-1 cells. Conversely, ΔNp63α knockdown had an opposite effect in T3M4 cells. ΔNp63α potentiated EGF-mediated activation of ERK, Akt and JNK signaling. Chromatin immunoprecipitation and functional reporter assays demonstrated that ΔNp63α activated EGFR transcription. 14-3-3σ transcription was also positively regulated by ΔNp63α and we have previously shown that 14-3-3σ contributes to chemoresistance in pancreatic cancer cell lines. Conversely, shRNA-mediated knockdown of 14-3-3σ led to abrogation of the ΔNp63α effects on cell proliferation and invasion. Thus, p53 homolog ΔNp63α enhances the oncogenic potential of pancreatic cancer cells through trans-activation of EGFR and 14-3-3σ.